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^J Abstract. We study the circular motion of charged test particles in the 

gravitational field of a charged mass described by the Reissner-Nordstrom 
spacetime, focusing on the physical differences between black holes and naked 
singularities. We perform the most complete classification of circular orbits for 
different sets of the main physical parameters, and study numerically the behavior 
of the angular momentum and energy of the charged test particle. Our analysis 
shows in an alternative manner that the behavior of circular orbits can be used 
^\ to distinguish between black holes and naked singularities. 

o\ 

<N 

• 

O 
CO 






> 



PACS numbers: 04.20.-q, 04.70.Bw, 04.40.Dg 

1. Introduction 

• • The Reissner-Nordstrom (RN) line element 

A 2 

ds 2 = -^dt 2 + r -dr 2 + r 2 (d9 2 + sin 2 Odcf) , (1) 

in standard spherical coordinates, describes the background of a static gravitational 
source of mass M and charge Q, where A = (r — r+)(r — r-) and r± = M± \/M 2 — Q 2 
are the radii of the outer and inner horizon, respectively. 

The motion of a test particle of charge q and mass \i moving in a RN background 
(HI) is described by the following Lagrangian density: 

C=±g a px 0l xP + eA a x 01 , (2) 

where A a are the components of the electromagnetic 4-potential, the dot represents 
differentiation with respect to the proper time, and the parameter e = q/ja is the 
specific charge of the test particle. 
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Since the Lagrangian density ([2| does not depend explicitly on the variables t 
and (/>, the following two conserved quantities exist 

dC /A. eQ\ E 

at \ r r / M 

P<f> = ^y =^sm 2 6^= -, (4) 

a0 M 

where L and E 1 are respectively the angular momentum and energy of the particle as 
measured by an observer at rest at infinity. 

On the equatorial plane = 7r/2, the motion equations can be reduced to the 
form r 2 + V 2 = E 2 / fi 2 which describes the motion of a test particle inside an effective 
potential V. Then, it is convenient to define the potential 




V ± = — = ^±\ 1+-2-- I" — + ^ (5) 



which corresponds to the value of E/jjl at which the (radial) kinetic energy of the 
particle vanishes [TJ |2j |3j [4], i. e. it is the value at which r is a "turning point" 
(V = E/fi). The effective potential with positive (negative) sign corresponds to the 
solution with 

lim E± = ±/i. 

r— )>oo 

Notice that in general E+ > E_ and E+(L, e, r) = —E_(L, — e, r). 

In the limiting case of vanishing charge, the effective potential reduces to 




V = +\[i + -o^)[l- — + ^r) ■ (6) 



This case was analyzed previously in [5j [6]. In this work, we first found that the 
stability properties of neutral test particles strongly depend on the nature of the 
central source. Indeed, in the case of a RN black hole there exists a minimum radius 
at which the orbit is stable, and outside this radius all the orbits are stable so that 
there exists only one region of stability. In the case of a naked singularity, the situation 
is completely different; the region of stability splits into two nonconnected regions so 
that a zone appears inside which no stable circular orbit can exist. This means that 
the stability properties of circular orbits could, in principle, be used to differentiate 
between a black hole and a naked singularity. 

In the general case of charged test particles, the presence of the additional 
term ^ changes completely the behavior of circular orbits, and leads to several 
possibilities which must be analyzed separately for black holes and naked singularities. 
In particular, it turns out that even in the case of the positive solution E+ and eQ < 0, 
negative energy states can exist [7] and also [8j [9j [TUl [TTJ [I2j [I3j [14]. In [15], we 
studied the spatial regions where circular motion is allowed around either black holes 
or naked singularities. As an example of the action of repulsive gravity, we found all 
the circles at which a particle can have vanishing angular momentum. We showed 
that it is possible to use the geometric structure of stable accretion disks, made of 
only test particles moving along circular orbits around the central body, to clearly 
distinguish between black holes and naked singularities (see also [16j[T7] and [18] [19]). 
For a recent discussion on the matter of naked singularity and black hole physics see 
[20] [21] and [22] [23] [24] . Also, there is extensive literature about particle motion in 
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the RN background. We refer in particular on [TS] and reference therein. More recent 
generalizations of previous analysis are given in [25 and [26 . 

In this work, we deepen the analysis performed in [15] in order to obtain the 
complete classification of circular motion around a RN black hole in section [2] and 
around a RN naked singularity with < e < 1, in section pi~T| and with — 1 < e < 0, 
in section 4.2 Conclusions close the paper in section ([5| 



2. Circular motion 

The effective potential (|5| regulates the circular motion of charged test particles. Only 
the special case of the positive solution V+ will be considered here in order to be able to 
compare our results with those obtained in the case of neutral test particles analyzed 
in [5j [6] . The extrema of the function V+ defined by the relations 

determine the radius of circular orbits and the corresponding values of the energy E 
and the angular momentum L. In the following analysis, we drop the subindex (+). 
Solving (|7|) with respect to L, we find the specific angular momentum 

(L ± ) 2 _ r 2 
/i 2 2S 2 

of the test particle on a circular orbit of radius r, where 

E = r 2 - 3Mr + 2Q 2 , (9) 

and the corresponding energy: 



2(Mr - Q 2 )S + e 2 Q 2 A ± QA ^e 2 (4S + e 2 Q 2 )l , (8) 



E± eQ A^2£ + e 2 Q 2 ±QV?(4^ 



• (io) 

/i r V2r\Z\ V J 

An interesting particular orbit is the one for which the particle is located at rest as seen 
by an observer at infinity, i.e., with L = 0. These "orbits" are therefore characterized 
by the following conditions 

dV 

L = 0, - = 0. (11) 



[27] . Solving (11) for Q ^ and e/0, we find the following radius 




_ (e 2 - 1) Q>M , , /v „ _ 

s ~ e 2 Q 2 -M* 1/ ' — — 2 , l j 

which is a generalization of the classical radius r* = Q 2 jM of a mass M with charge 
Q. A detailed analysis of the behavior of this radius was presented in [T5] . 

3. Circular motion around a RN black hole 

In this section, we present a classification of the circular orbits around a RN black 
hole (Q < M) by using certain values of the parameters L and e which follow from 
the conditions of circular motion. It can be considered as an alternative classification 
to that presented in [15 , in which the value of the ratio Q/M plays a central role. 
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20 r 




Figure 1. The charge parameters e\ (dashed curve), e (black solid curve), e n 
(dotdashed curve), and M/Q (gray curve) as functions of the charge-to-mass 
ratio of the RN black hole. 




Figure 2. The charge parameters Q n (gray curve), Q (black solid curve), and 
M/e (dashed curve) as functions of the charge-to-mass ratio e of the test particle. 



Here we are interested in the region outside the outer horizon r+. Let us introduce 
the parameters 



—L^5M 2 - 4Q 2 + V^M 2 - 24Q 2 , 



Cn = 



/3M 2 v/(9M 2 - 8Q 2 )M 2 



2 + e^ 1 + e z 

the behavior of which is depicted in figures [T] and [2] 
Moreover, let us introduce the radius 
, 3M . 1 



2 ±-V9M2-8Q 2 , 



(13) 
(14) 
(15) 

(16) 



which represents the limiting radius at which neutral particles can be in circular motion 
around a RN naked singularity [5_, the angular momentum parameter and 

9M 7 /„„ . /^tt^ — r^7\ . 2Q 2 M 2 



L 



" ~ 2e 2 Q 2 



(3M + a/9M 2 - 8Q 2 ) 



(1 + e 2 ) 
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Table 1. Characteristics of the circular orbits for charged test particles with 
charge-to-mass ratio e < 1 in a RN black hole. 



0< e< 1 






€ < 


Region 


Momentum 


Region 


Momentum 



<r <rZ 



r > r2 



L ■ 

L - 

L - 

L - 



■■±L± 

±L n 



r > r2 



±L+ 



2e 2 



27M + 5a/9M 2 -8Q 2 
which satisfies the relationship (dV/dr){L, 



3 (ZM + a/9M 2 - 8Q 2 ) e 2 ] , (17) 



+^ - 



) = 0, and the radius 



3M 

~2~ 



i^x/QM 2 



8Q 2 - Q 2 



where 



lim r. 



(18) 
(19) 

(20) 



We consider separately the two cases e < 1 and e > 1 and investigate the value of 
the allowed orbit's radius and the corresponding angular momentum. The results are 
summarized for e < 1 in Table [TJ and for e > 1 in Tables [3] and [4] Moreover, the 
characteristic parameters for all the different cases are listed in Table [2] 

We can summarize the results as follows. For e < 1, Table [I] the infimum circular 
orbital radius is located at r+ (r+) for e < (0 < e < 1). The situation is much more 
complicated for the case (0 < e < 1). In this case, the electromagnetic interaction is 
repulsive, and is balanced by the attractive gravitational component. Different values 
of the angular momentum are possible in the set {±L±, ±£ n }- The region in which 
circular orbits are allowed is therefore split by the radius r+ and r^ . 

The case e > 1 should be considered apart. This corresponds to the real case 
of charged elementary particles, like electrons, protons and ions, orbiting around a 
charged source. In the case eQ > with e > 1, the circular motion dynamics is 
determined by the classification given in[3j for a fixed source charge-to-mass ratio, or 
in[4j for a fixed particle charge-to-mass ratio. 

First, we can recognize nine different situations listed in Table pi). There are 
orbits with angular momentum L = ±L±, L = ±L n or even L = 0. This last case 
occurs when the particle is located at rest with respect to an observer at infinity. 
Clearly, these particular "orbits" are the consequence of the full balance between 
the attraction and repulsion among test particles and sources. The infimum circular 

'£. Comparing with the case 



orbital radius can be r 



rt or also r 



e < 1, we can clearly see that the main difference between the two cases consists in 
the presence of a supremum circular orbital radius at r = r+ or at r = rf. This 
means that the repulsive electromagnetic effect, predominant at larger distances, does 
not allow circular orbits around the source. The classifications of Table |3] and Tabled] 
are alternative and equivalent. In the first one, if we fix the charge-to-mass ratio of 
the black hole, and we move towards increasing values of the particle charge-to-mass 
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Table 2. Values of the angular momentum L which are possible in different 
regions of the radial coordinate r. 



Case I 




Case II 




Case III 




Region 


Momentum 


Region 


Momentum 


Region 


Momentum 


r = rf 

r = rf 
(rf , oo) 


L = ±L_ 
L = ±L± 
L = ±L n 

L = ±L_ 


r = rf 
(r+,r+) 
r = rf 
(rf,rf) 


L = ±L- 
L = ±L± 
L = ±L n 

L = ±L- 


r = r~f~ 
r = Tj 


L = ±L- 
L = ±L± 
L = 


Case IV 




Case V 




Case VI 




Region 


Momentum 


Region 


Momentum 


Region 


Momentum 


r = rf 
(rf,rf) 

r = rf 
(r+,r+) 


L = ±L_ 

L = ±L± 

L = ±L + ; L = 

L = ±L + 


r = rf(rf) 
(r z + ,r+)((r+,r+)) 


L = 

L = ±L + 


r = Tj 
(r+,r+) 

r = r^ 
(r|,r+) 
r = rf 


L = ±L- 
L = ±L± 
L = ±L n 
L = ±L- 
L = 


Case VII 




Case VIII 




Case IX 




Region 


Momentum 


Region 


Momentum 


Region 


Momentum 



r = rf(r£) 


±L- 


r = r + 


L = 


r = rf 


L = 


(r+,r+) 


L = ±L± 


(r+,r+) 


L = ±L + 


(r+,r+) 


L = ±L + 


r = rf 


L = 


r = r^j 


L = ±L n 







Table 3. Classification of the circular orbits of a charged test particle with 
charge-to-mass ratio e > 1 in terms of the value of the ratio Q/M of a RN black 
hole. 



Q/M e (0, 1/2] 


Q/M e [1/2, x/13/5) 


q/m e [Vis/6, VVS) 


Q/M e [v^Ts, i) 


Region 


Case 


Region 


Case 


Region 


Case 


Region 


Case 


1< e<2 


I 


1< e < M/Q 


I 


1< e < M/Q 


I 


1< e < M/Q 


/ 


2 < e < M/Q 


I 


M/Q < e < 2 


II 


M/Q <e<e n 


II 


M/Q <6<6 n 


II 


M/Q <c<c n 


VI 


2 < e < e n 


VI 


e = e n 


III 


e = e n 


III 


e = e n 


VII 


e = e n 


VII 


e n < e < 2 


IV 


e n < e < e 


IV 


e n < e < e 


IV 


e n < e < e 


IV 


2 < e < e 


IV 


e < e< 2 


V 


€ = € 


VIII 


e = e 


VIII 


e = e 


VIII 


e>2 


IX 


e > I 


IX 


e > e 


IX 


e > I 


IX 







ratio, we find four classes of objects characterized by 



I 

II 
III 
IV 



Q/M e (0,1/2], 
Q/M e (1/2, VI3/5] , 
Q/M e (\/l3/5, a/273) , 
pe]#,l), 



(21) 
(22) 
(23) 
(24) 



respectively. Thus, for a fixed charge of the spacetime, and following Table [3] and 
Table ([2|, we can trace a picture of the dynamical properties of that spacetime. An 
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Table 4. Classification in terms of the charge-to-mass ratio e > 1 of the circular 
orbits of a charged test particle moving in the field of a RN black hole with mass 
M and charge Q. 



1< € < 2 




e> 2 




Region 


Case 


Region 


Case 


< Q < M/e 


I 


< Q < M/e 


/ 


M/e <Q <Q n 


II 


M/e <Q<Q n 


VI 


Q = Qn 


III 


Q = Qn 


VII 


Qn <Q <Q 


IV 


Qn <Q <Q 


IV 


Q<Q<M 


V 


Q = Q 


VIII 






Q <Q <M 


IX 



Q=0.1M 



3.2 
3.0 
2.8 
| 2.6 
2.4 
2.2 
2.0 



Q=0.1M 



Ml 



20 40 




Figure 3. Black hole case with Q = 0.1M. Left panel: The radiai rj, r^~ and 
r~/ are plotted as functions of the test particle charge-to-mass ratio e. Shaded 
regions are forbidden (Right panel). 



example is provided in figure [3] and figure [4j where the case of a spacetime of charge- 
to-mass ratio Q = 0.1M is illustrated. On the other hand, in Table |1J if we fix the 
charge-to-mass ratio of the test particle, and move towards increasing values of the 
charge-to-mass ratio of the source, we find only two different cases: e < 2 and e > 2. 
Using this classification, we can determine the orbital radius followed by the selected 
particle in the fixed spacetime. 



4. Circular motion around a RN naked singularity 

Equations u?\ govern the circular motion around a RN naked singularity (Q < M). 
Because (pi) and (10) define the angular momentum L± and the energy E± in terms 



of r/M, Q/M, and e, it is necessary to investigate several intervals of values. To this 
end, as in [15J, it is useful to introduce the following notation 



Q 



_ V9M 2 



e± = 



V2Q 



5M 2 ± 4Q 2 + V25M 2 - 24Q 2 



(25) 
(26) 
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L + 
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1 


2.0 
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Figure 4. Black hole case with Q 

radiai rf , 

ratio e in the region [—10,30]. 



0.1M. Shaded regions are forbidden. The 



e± = 



V^Q 



3M 2 - 2Q 2 =b M^9M 2 



(27) 



where e_ = e and e + = e n as defined in (13) and (14), respectively. For completeness, 
and following [15 , we reproduce here in figure plthe behavior of these parameters in 
terms of the ratio Q/M > 1. It follows from figure [5] that it is necessary to consider 
the following four intervals: 

I : Q/M €(1, 5/(2>/6)], (28) 

II : Q/M e (5/(2>/6), (3\/6)/7] , (29) 

III : Q/M e ((3>/6)/7, v/9/8] , (30) 

IV: Q/M g [V9/8,oo). (31) 

This analysis completes the classification briefly presented in [15], and now contains 
all the different cases that are allowed by the values of e. 



For e > 0, in general, circular orbits exist in the region r > 



= n2 



Q 2 /M. In 



e > 1 and M < Q ■ 
- L+ in the interval r~ 



C \J 9/8M circular orbits exist with angular 
< r < r+ while for Q > y/9/SM no circular 



particular, for 

momentum L - 

orbits exist [15]. The energy and angular momentum of circular orbits diverge as 

approaches the limiting orbits at r^: the boundary r = r+ in this case corresponds 

to a lightlike hypersurface. Finally, the energy of the states is always positive. 

Since for eQ > the Coulomb interaction is repulsive, the situation characterized 
by the values for Q > y/9/8M and e > 1 corresponds to a repulsive electromagnetic 
effect that cannot be balanced by the attractive gravitational interaction. 
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" 1.00 1.01 1.02 1.03 1.04 1.05 1.06 

Q/M 

Figure 5. The charge parameters ei (black solid curve), e_ (gray solid curve), 
e+ (dashed curve), e_ (dotted curve), and e+ (dotdashed curve) as functions 
of the charge-to-mass ratio of the RN naked singularity. The special lines 
Q/M = 5/ (2^) « 1.02, Q/M = ?>Vq/7 « 1.05, and Q/M = yW^ « 1-06 
are also plotted. 



For e < —1, the contribution of the electromagnetic interaction is always 
attractive. Hence, the repulsive force necessary to balance the attractive effects of 
the gravitational and Coulomb interactions can be generated only by the RN naked 
singularity. In particular, for e < —1 and for Q > y^9/SM circular orbits with L = L+ 
always exist for r > 0. For M < Q < ^9/8M, circular orbits exist with L = L+ in 
< r < r ~ and r > r+ . Charged test particles with e < — 1 can move along circular 
orbits also in the region (0,r*]. The value of the energy on circular orbits increases 
as r approaches r = 0, and the angular momentum, as seen by an observer located at 
infinity, decreases as the radius of the orbit decreases. In the region M < Q < y/9/8M, 
two limiting orbits appear at r^ (similar to the neutral particle case [5]). 



4-1. Circular motion of particles with < e < 1 in a RN naked singularity 

In this case, it is necessary to consider separately all the possible values of the charge 
parameters e/, e± and e± in all the regions given in (28JJ31) which, according to figure 
[5j follow for the values of the ratio Q/M. 

In the first region, 1 < Q/M < 5/(2\/6), we analyze the energy and angular 
momentum of charged test particles for the following cases: 

(i) M <Q < 5/(2 V6)M and < e < 1_: See figure [HJ 

(ii) M <Q< 5/(2 V6)M and !_ < e < e_: See figure [7 

(hi) M <Q< 5/(2 V6)M and e_ < e < e + : See figure | 

(iv) M <Q < 5/(2 VE)M and e + < e < q: See figure [9] 

(v) M <Q < 5/(2 >/6)M and e z < e < e + : See figure [To} 

(vi) M <Q< 5/(2 >/6)M and ^ + < e < M/Q: See figure [TT] 

(vii) M <Q< 5/(2 v / 6)M and M/Q < e < 1: See figure [l§ 

In the second interval, 5/(2\/6) < Q/M < (3\/6)/7, we analyze the energy and 
angular momentum of charged test particles for the following cases: 



(i) 5/(2 VE)M < Q < (3V6/7)M and < e < !_: See figure 13 
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1.040 
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1.041 1.042 




1.957 1.958 




Figure 6. Case: M < Q < 5/(2v / 6)M and < e < e_. Parameter choice: 



1.01M and e = 0.12. Then e_ = 0.14639, rj = 1.03789M, r~ 



1.04196M, 



r+ = 1.95804M, r~ = 1.04599M, and r+ = 1.95401M. Circular orbits exist with 
angular momentum L = L~ (gray curve) and energy E = E~ (black curve) in 



< r < r^ 



at r ■ 



(upper left plot); L = L± in r 7 < r < r l (upper 



right plot) and in r^<r< rZJT (bottom left plot); L = L~ in r > rZJT (bottom 
right plot). The angular momentum L+ is represented by a gray dotted curve 
and the energy E+ by a black dashed curve. 



(ii) 5/(2 VE)M < Q < (3V6/7)M and e_ < e < e t : See figure 

(iii) 5/(2 v / 6)M < Q < (3V®/7)M and e t < e < e + : See figure 

(iv) 5/(2V6)M < Q < (3v / 6/7)M and 1+ < e < M/Q: See figure 

(v) 5/(2V / 6)M < Q < (SV^/TjM and M/Q < e < 1: See figure 



16 



17 



In the third interval, 3^/7 < Q/M < \/9/8, we analyze the energy and angular 
momentum of charged test particles for the following cases: 

(i) (3\/6/7)M< 



18 



19 



< V9/8M and < e < e^: See figure 
(ii) (3 v / 6/7)M < Q < y/9/SM and e/ < e < 1_: See figure 

(iii) (3\/6/7)M < Q < v / 978 M and l - < e < ^+ : See fi g ure 
(iv) (3v / 6/7)M < Q < y/9/SM and 1+ < e < M/Q: See figure 
(v) (3\/6/7)M < Q < V^M and M/Q < e < 1: See figure 



20 



21 



22 



In the region Q > y/9/SM^ the only case in which circular orbits exist is for 
< e < M/Q. See figure [23] In general we can summarize the situation as follows. 

There is always a minimum radius r m i n at which circular motion is allowed. At 
the radius r~ the energy of the test particle diverges, indicating that the hypersurface 
r = r~ is lightlike. In the simplest case, there is a minimum radius r m ^ n so that 
circular orbits are allowed in the infinite interval ]r m i n , oo). Otherwise, this region is 
split by a lightlike hypersurface situated at r+ > r min . 
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1.0505 1.0510 1.0515 1.0520 1.0525 1.0530 
r/M 



1.948 1.950 1.952 1.954 1.956 1.958 
r/M 




Q=1.01M 
e=0.2 



6 

r/M 



Figure 7. Case: M < Q < 5/(2v / 6)M and e_ < e < e_. Parameter choice: 
Q = 1.01M and e = 0.2. Then e_ = 0.14639, e_ = 0.313219, rf = 1.0505M, 
r~ = 1.04196M, rf = 1.95804M, r~ = 1.05323M, and rf = 1.94677M. Circular 
orbits exist with angular momentum L = L+ (gray curve) and energy E = E+ 
(black curve) in rvJT < r < rj (upper left plot); L = L± in r J < r < r^~ (upper 
right plot) and r^ + < r < rf (bottom left plot); L = L~ in r > rf (bottom right 
plot). The angular momentum L~ (gray dotted curve) and the energy E~ (black 
dashed curve) are also plotted. 




1.90 1.91 1.92 1.93 1.94 1.95 
r/M 



Figure 8. Case: M < Q < 5/(2 v / 6)M and e_ < 
Q = 1.01M and e = 0.5. Then e_ = 0.313219, e+ = 



: < e+. Parameter choice: 

0.896311 rf = 1.11104M, 

1.04196M, r+ = 1.95804M, r," = 1.11784M and r+ = 1.88216M. Circular 



orbits exist with angular momentum L = L+ (gray curve) and energy E = E. 



(black curve) in r 7 < r < rf (left plot); L = at r : 
(center plot); L = L~ in r > rf (right plot). 



L± in rf < r < n 



:t 
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Figure 9. Case: M < Q < 5/(2v / 6)M and e + < e < e t . Parameter choice: 



1.01M and e = 0.902. Then e t = 0.907, e+ 



0.8963, r] 



1.04196M, r+ = 1.95804M, 



1.45192M, and r+ 



= 1.44942M, 
1.548077M. 

Circular orbits exist with angular momentum L = L+ (gray curves) and energy 
E = E+ (black curves) in r^ < r < rf~ (upper left plot); L = L± in r j" < r < rj~ 
(upper right plot) and r^~ < r < ?%Jr (bottom left plot); L = L~ in r > rZjr (bottom 
right plot). 





Figure 10. Case: M < Q < 5/(2v / 6)M and e L < e < e+. Parameter choice: 
Q = 1.01M and e = 0.92. Then ei = 0.907, e+ = 0.958884, r+ = 1.528942M, 
r~ = 1.04196M, and r^T = 1.95804M. Circular orbits exist with angular 
momentum L = L+ (gray curves) and energy E = E+ (black curves) in 
rj < r < r]~ (left plot); L = L± in rj~ < r < rZJT (center plot); and L = L~ in 
r > rZJT (right plot). 
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Figure 11. Case: M < Q < 5/(2 v / 6)M and e+ < e < M/Q. Parameter choice: 

Q = 1.01M and e = 0.98. Then 1+ = 0.958884, r+ = 3.37999M, r~ = 1.04196M, 
and rj = 1.95804M. Circular orbits exist with angular momentum L = L+ (gray 
curves) and energy E = E+ (black curves) in r^ < r < rZJT (left plot) and L = L~ 
in r > rj (right plot). 




1.2 



1.4 1.6 

r/M 



Figure 12. Case: M < Q < 5/(2 v / 6)M and M/Q < e < 1. Parameter 
choice: Q = 1.01M and e = 0.995. Then e + = 0.958884 r+ = 0.422853M, 
r~ = 1.04196M, and r^ = 1.95804M. Circular orbits exist with angular 
momentum / = L+ (gray curve) and energy E = E+ (black curve) in r^ < r < rZJT . 



Another possible structure is that of a spatial configuration formed by two 
separated regions in which circular motion is allowed in a finite region filled with 
charged particles within the spatial interval (r min = r~ ^r^ax — r+). This region is 
usually surrounded by an empty finite region in which no motion is allowed. Outside 
the empty region, we find a zone of allowed circular motion in which either only neutral 
particles or neutral and charged particles can exist in circular motion. 
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1.10 1.11 1.12 1.13 



1.142 1.143 1.144 1.145 1.146 1.147 




Figure 13. Case: 5/(2 v / 6)M < Q < (3V6/7)M and < e < e_. Parameter 
choice: Q = 1.03M and e = 0.12. Then e_ = 0.2772, r+ = 1.0935M, 
r~ = 1.1495M, r+ = 1.858M, r~ = 1.14732M, and r+ = 1.8526804M. Circular 
orbits exist with angular momentum L = L~ (gray curves) and energy E = E~ 
(black curves) in r J < r < r^ (upper left plot); L = at r = rj; L = L± in 
Tj<r< r~^~ (upper right plot) and ry < r < r^ (bottom left plot); L = L~ in 
r > rZJT (bottom right plot). 




Figure 14. Case: 5/(2y/6)M < Q < (3 v / 6/7)M and e_ < e < q. Parameter 
choice: Q = 1.03M and e = 0.5. Then £_ = 0.2772, e t = 0.695243, r+ = 1.237M, 
r~ = 1.1495M, r+ = 1.858M, r~ = 1.25122M, and rf = 1.74878M. Circular 
orbits exist with angular momentum L = L+ (gray curve) and energy E = E+ 
(black curve) in r^ < r < rs (upper left plot); L = at r = rj; L = L± in 
rj<r< r~^ (upper right plot) and ry < r < r^ (bottom left plot); L = L~ in 
r > rZJT (bottom right plot). (Reproduced from |15|.) 
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Figure 15. Case: 5/(2 v / 6)M < Q < (3 v / 6/7)M and e t < eQ < e+. Parameter 
choice: Q = 1.03M and e = 0.8. Then e+ = 0.866828 and e t = 0.695243, 
r+ = 1.58116M, r~ = 1.1495M, and r+ = 1.858M. Circular orbits exist with 
angular momentum L = L+ (gray curve) and energy E = E+ (black curve) in 
r^ < r < rf (left plot); I = L± in r t < r < r^ (center plot); and L = L~ in 




Figure 16. Case: 5/(2 v / 6)M < Q < (3v / 6/7)M and e + < e < M/Q. Parameter 
choice: Q = 1.03M and e = 0.92. Then 1+ = 0.866828, r+ = 2.52M, 
r^ = 1.1495M, and r^j" = 1.858M. Circular orbits exist with angular momentum 
L = L+ (gray curve) and energy E = E+ (black curve) in r^ < r < rZJT ( left 
plot) and L = L~ in r > rj~; L = at r = rj~ (right plot). 




Figure 17. Case: 5/(2 v / 6)M < Q < (3 v / 6/7)M and M/Q < eQ < 1. 

Parameter choice: Q = 1.03M and e = 0.92. Then e + = 0.866828, r+ = 2.52M, 

r^ = 1.1495M, and r^ = 1.858M. Circular orbits exist with L = L+ (gray 
curve) and E = E+ (black curve) in r^ < r < rZJT. 
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Q=1.05M 
e=0.12 




1.16 1.18 1.20 1.22 1.24 1.26 1.28 
r/M 




1.290 1.292 1.294 1.296 

r/M 




1.704 1.706 1.708 

r/M 



Figure 18. Case: (3y/6/7)M <Q< ^/9/8M and < e < e t . Parameter choice: 
Q = 1.05M and e = 0.12. Then e t = 0.404, r+ = 1.1468M, r~ = 1.28787M, 
r+ = 1.71213M r~ = 1.29744M, and r+ = 1.70256M. Circular orbits exist with 
angular momentum L = L~ (gray curve) and energy E = E~ (black curve) in 



< r < r. 



■■ at r - 



(upper left plot); L = L± in r^ < r < r l (upper 



right plot) and r^~ < r < r^ (bottom left plot); L = L 
plot). 



in r > rj (bottom right 




1.2850 1.2855 1.2860 1.2865 1.2870 1.2875 
r/M 




1.3 1.4 



1.6 1.7 




Figure 19. Case: (3v / 6/7)M < Q < ^/9/8M and e L < e < e_. Parameter choice: 
Q = 1.05M and e = 0.405. Then e t = 0.404, e_ = 0.41, r+ = 1.28M, r~ = 



1.28787M, and r^ = 1.71213M. Circular orbits exist with angular momentum 

L = L~ (gray curve) and energy E = E~ (black curve) in r t < r < r^T; L = at 

+ 

7 



r = rj (left plot); L = L± in r 7 
(right plot). 



< r < r^ (center plot); and L = L in r > r. 
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1.35 1.40 



1.55 1.60 1.65 1.70 



Figure 20. Case: (3y/6/7)M < Q < ^/9/8M and e_ < e < e+. Parameter 
choice: Q = 1.05M and e = 0.5. Then 1+ = 0.74, e_ = 0.41, rf = 1.35M, r~ = 
1.28787M, and rZJT = 1.71213M. Circular orbits exist with angular momentum 
L = L+ (gray curve) and energy E = E+ (black curve) in r^ < r < rf (left 
plot); L = L± in r f < r < rZJT (center plot); and L = L~ in r > r^~ (right plot). 




Figure 21. Case: 
Q = 1.05M and e = 



(3^6/7)M < 
: 0.92. Then e+ 



/978M and 1+ < e < M/Q. choice: 



4.437M, r" 



1.28787M, and 



r^ = 1.71213M. Circular orbits exist with angular momentum L = L+ (gray 
curve) and energy E = E+ (black curve) in r^j < r < r^ (left plot) and L = L~ 



in r > rf; L = at r 



(right plot). 




Figure 22. Case: (3 v / 6/7)M < Q < ^/9/8M and M/Q < e < 1. 



choice: Q = 1.05M and e = 0.92. Then r~ 



1.28787M and r+ 



Parameter 
1.71213M. 



Circular orbits exist with angular momentum L = L+ (gray curve) and energy 
E = E+ (black curve) in r^ < r < rZJT. 
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9.0 9.5 10.0 



Figure 23. Case:0 < e < M/Q and Q > ^/9/8M. Parameter choice: Q = 2M 
and e = 0.2. Then rj = 6.18767M. Circular orbits exist with angular momentum 
L = L~ (gray curve) and energy E = E~ (black curve) in r > r t and L = at 
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Figure 24. Case: Q > ^/9/8M and -1 < e < —M/Q. Parameter choice: 
Q = 2M and e = -0.8. Then r+ = 1.2086M. Circular orbits exist with angular 
momentum L = L+ (gray curve) and energy E = E+ (black curve) in r > rj. 
For r = r+, L = 0. 




Figure 25. Case: Q > ^/9/8M and e = -M/Q. Parameter choice: Q = 2M 
and e = —0.5. Then rj = 1.2086M. Circular orbits exist with angular 
momentum L = L+ (gray curve) and energy E = E+ (black curve) in the region 
r > Q 2 /(2M). For r = Q 2 /(2M), L = 0. 



4-2. Circular motion of particles with —l<e<0ina RN naked singularity 

For this range of values of the ratio e, we analyze the angular momentum and the 
energy in the following cases: 



> V9/8M and -1 < e < -M/Q: See figure 

> V^M and e = -M/Q: See figure 



24 



25 



(i) 

(ii) 

(hi) Q > y/9/SM and -M/Q < e < 0: See figure 

(iv) M <Q < V/9/8M and -1< e < -M/Q: See figure 

(v) M <Q< V^M and e = -M/Q: See figure 



26 



27 



28 



(vi) M <Q< V9/8M and -M/Q < e < 0: See figure 



29 



In general two different configurations are allowed. For Q > y/9/SM a continuous 
region appears from a minimum radius r m i n G {r^, Q 2 /(2M)} to infinity in which 
circular orbits are allowed. 
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Figure 26. Case: Q > ^/9/8M and -M/Q < e < 0. Parameter choice: Q = 2M 
and e = —0.2. Then r~^ = 2.9M. Circular orbits exist with angular momentum 
L = L+ (gray curve) and energy E = E+ (black curve) in r > rj . For r = 7\T , 
L = 0. 




Figure 27. Case: M < Q < ^/9/8M and -1 < e < -M/Q. Parameter 
choice: Q = 1.05M and e = 0.98. Then r+ = 0.427878M, r~ = 1.28787M, and 
r^j = 1.71213M. Circular orbits exist with angular momentum L = L+ (gray 
curve) and energy E = E+ (black curve) in r J < r < r^ (left plot) and in r > r^ 



for (right plot). For r - 



r> + 



L = 0. (Reproduced from |15|.) 



For M < Q < y/9/SM there is a non connected region (r m ^ n ,oo) inside which 
there is a forbidden region (r~ ,r+). The configuration is therefore composed of two 
disconnected regions. 



5. Conclusions 



In this work, we explored the motion of charged test particles along circular orbits 
in the spacetime described by the Reissner-Nordstrom metric. A detailed discussion 
of the dynamics for the black hole and naked singularity cases has been performed. 
Circular orbits have been classified in detail for a complete set of cases. 

We adopt the effective potential approach to study test particle orbits in the 
Reissner-Nordstrom spacetime, focusing on the equatorial orbits. We explore the 
morphology of the orbital regions. This analysis leads to a clear differentiation between 
naked singularities and black holes. A remarkable implication of this description is 
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6 
r/M 



Q=1.05M 
e=-l/(1.05) 



Figure 28. Case: M < 
Q = 1.05M. Then r~ 



/9/8M and e 



-M/Q. Parameter choice: 



1.28787M and r+ = 1.71213M. Circular orbits exist 
with angular momentum L = L+ (gray curve) and energy E = E+ (black curve) 
in Q 2 /(2M) < r < r~ (left plot) and in r > r+ (right plot). For r = Q 2 /(2M), 
L = 0. 




Figure 29. 

choice: Q - 



Case: M < 
1.05M and e 



< y^M and -M/Q < e < 0. 



-0.2. Then rZ 



1.28787M, 



Parameter 
1.71213M, 



and r s = 1.03487M. Circular orbits exist with angular momentum L = L+ 
(gray curve) and energy E = E+ (black curve) in rj < r < r^ (left plot) and in 
r > ?%jr (right plot). For r = r^~ , L = 0. 



that the circular orbit configuration in the black hole case is not allowed in the naked 
singularity regime. Instead, the study of the circular orbits trace out a possible way 
to distinguish the two physical situations. 

In a series of previous works [5j |6j [15], we analyzed the dynamics of the RN 
spacetime, and studied the motion of neutral and charged test particles, by using the 
effective potential approach. We showed that in the case of charged test particles the 
term eQ/r drastically changes the behavior of the effective potential. The study shows 
the existence of stability regions whose geometric structure clearly distinguishes naked 
singularities from black holes (see also [16j[l7] and p~8j[19]). In [15], in particular, we 
studied the spatial regions of the RN spacetime where circular motion is allowed 
around either black holes or naked singularities. We showed that the geometric 
structure of stable accretion disks allows us to clearly distinguish between black holes 
and naked singularities. In this work, we presented the complete classification of 
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circular motion around a RN black hole, and around a RN naked singularity with 
< e < 1 and -1 < e < 0. 

Clearly, this analysis could be used to construct an accretion disk with 
disconnected rings made of test particles. A precise characterization of matter 
configurations surrounding a charged compact astrophysical object can account for 
significant astrophysical processes observed in the electromagnetic band, like the 
jet emissions. Therefore understanding the dynamics around compact objects is 
important for the understanding of the accretion disk phenomena and the classification 
of their general properties. 

In the naked singularity case, gravity can assume a repulsive nature that produces 
a complicate picture of dynamics around the source. We distinguish four regions of 
charge-to-mass ratio values of the source that characterized these configurations. In 
each region, different cases can occur: an infinite, continuum, region or otherwise a 
disconnected region. In the black hole case the circular orbits configuration strongly 
varies if e < 1 or e > 1. An infimum radius always appears, and in the case e > 1 a 
supremum circular orbits radius can appear. 
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